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SUMMARY 


An investigation was made to determine the accuracy with which the 
lateral flight motions of a swept-wing airplane could he predicted from . 
experimental stability derivatives determined in the 6 -foot -diameter 
rolling-flow test section and 6- by 6-foot curved-flow test section of 
the Langley stability tunnel. In addition determination of the significance 
of including the nonlinear aerodynamic effects of sideslip in the 
calculations of the motions was desired. All experimental aerodynamic 
data necessary for prediction of the lateral flight motions are presented 
along with a number of comparisons between flight and calculated motions 
caused by rudder and aileron disturbances. 

In general, the agreement between the calculated and measured motions 
of the airplane considered was good when the effects of all control 
movements were taken into account. The greatest disagreement occurred 
at lift coefficients where Reynolds number effects on the experimental 
derivatives would be expected to be high, which for the case considered 
was for lift coefficients above about 0.8 when wing slots were used. The 
nonlinear effects of sideslip for this airplane . were not very significant 
for the motions considered, which generally involved sideslip angles less 
than 10°. 


INTRODUCTION 


For the past few years, numerous investigations have been made in 
the Langley stability tunnel to determine the effects of geometric 
variables on the static-, rolling-, yawing-, and pitching-stability 
derivatives of various airplane configurations. (See references 1 to 4.) 
In the past, however, none of the experimental data has been compared 
with data obtained in flight to determine its relative worth. The 
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purpose of the present paper. Is to determine the applicability of the 
experimental stability derivatives to the prediction of the lateral 
disturbed motions of an airplane in flight. The equations used for 
calculating the motions are given in the appendix. 

A ^-scale model of a swept-wing version of a conventional fighter 
airplane, which was selected because of the large amount of flight data 
available (see references 5 and 6), was tested in the 6-foot-diameter 
rolling-flow test section and 6- by 6-foot curved-flow test section of 
the Langley stability tunnel to determine all the stability derivatives 
which are usually considered necessary to calculate the lateral motions 
arising from a disturbance caused by the rudder or the ailerons. 
Comparisons have been made between the flight and calculated lateral 
motions for a wide range of conditions in gliding flight. 

A few calculations have been made to determine the effects of 
nonlinear variations of the aerodynamic forces and moments with the 
angle of sideslip. 


COEFFICIENTS AND SIMBOIS 


The results of the wind-tunnel tests are presented as standard NACA 
coefficients of forces and moments. Moment coefficients are referred 
to a center of gravity located at 21.8 percent of the mean aerodynamic 
chord. The wind-tunnel data and motion calculations are referred to 
the stability axes, which are a system of axes having their origin at 
the center of gravity, and in which the Z-axis is in the plane* of symmetry, 
and perpendicular to the relative wind, the X-axis is in the plane of 
symmetry and perpendicular to the Z-axis, and the Y-axis is perpendicular 
to the plane of symmetry. The positive directions of the stability axes 
and of angular displacements of the airplane and control surfaces are 
shown in figure 1. 

The coefficients and symbols are defined as follows: 

C L lift coefficient (Lift/qS) 

Ct maximum lift coefficient 

-Toax 

C x longitudinal-force coefficient (x/qS) 

Cy lateral -force coefficient (Y/qS) 

Ci rolling -moment coefficient (L/qSb) 

C n yawing-moment coefficient (N/qSb) 
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moment of inertia about normal principal axis 


longitudinal force along X-axis 
lateral force along Y-axis 
normal force along Z-axis (Lift = -Z) 
rolling -moment about X-axis 
yawing -moment about Z-axis 
wing-tip helix angle, radians 
yawing-velocity parameter, radian measure 
rolling angular velocity about X-axis 
yawing angular velocity about Z-axis 
linear velocity of airplane along Y-axis 
free-stream velocity along X-axis 

calibrated airspeed; based on sea-level density of air 
angle of sideslip; P = -\|r in wind-tunnel tests / tan - ^- 


angle of attack' of root chord line 
angle of attack of thrust line - 1.2°) 
angle of yaw, degrees 

angle of incidence of stabilizer with respect to thrust line, 
positive when trailing edge is down 

control -surface deflection; measured in a plane perpendicular 
to hinge axis 


angle of sweepback, degrees 
free-stream dynamic pressure 


(t'j 


S 


wing area 
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wing span. 

aspect ratio ("b 2 /S) 
mass density of air 
time 

time to damp to half-amplitude 
period 
tail length 
Subscripts : 
a aileron 

r rudder 

f flap 

v vertical tail 

WMD-TUMEL TESTS 
Appratus and Model 


b 

A 

P 

T 

T l/2 

P 

l 


The experimental static-lateral-stability derivatives, rolling- 
stability derivatives, and yawing -stability derivatives were determined 
from tests conducted in the Langley stability tunnel in which roll ing 
or curved flight is simulated by rolling or curving the air stream about 
a rigidly mounted model. 

The tests were made on a conventional six -component balance system 
with the model mounted at the flight center of gravity which is at 
21.8 percent of the mean aerodynamic chord of the wing. 

The full-scale airplane (a swept-wing version of a conventional 
fighter) had the quarter-chord line of the wing, Just outboard of the 
intake ducts, swept back 35°. Some, of the pertinent airplane character- 
istics are given in table I. More details of the airplane may be 
obtained from references 5 and 6. 
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The —-scale model shown in figure 2 and in the photographs of 


figures 3 to 6 was constructed of laminated mahogany, finished in clear 
laquer, and ali surfaces were’ highly polished. The model propeller 
had three metal blades set at an angle of 28° at the 0.75 radius. All 
propeller-on tests were made with windmilling propeller. The model 
wing had a removable leading edge so that slats of 0 percent, kO percent 
and 80 percent of the swept span could be used interchangeably to 
simulate those of the full-scale airplane. The top surfaces of the 
slats were cast to the contour of the airfoil and the slats were 
extended by means of metal brackets which also act as fences to reduce 
spanwise flow along the slot. A cross section through the slot and 
slat is shown in reference 5. 




The wing had a plain trailing-edge flap with a chord of 15.1 percent 
of the wing chord measured perpendicular to the hinge axis. The gap was 
sealed for all tests. As in the case of the airplane the main wheels 
of the model were fixed for all tests; whereas the nose wheel and 
nose-wheel door3 were removed for all flaps -up tests. 


Shown in figure 2 are the two ventral fins tested on the model. 
The large ventral fin was used for all tests except a few with the 
80-percent-span slot configuration for which the small ventral fin 
was used. 


Tests and Test Conditions 

Trim tests .- Model trim lift coefficients of 0.33, 0.55, O.76, 
and 0.95 were selected as representative of those obtained in flight 
tests . The angle 'of incidence of the horizontal tail was measured 
with respect to the thrust line. 

In order to determine the trim angles of the horizontal tail for 
the previously mentioned lift coefficients, tests were made through the 
angle -of -attack range with the horizontal tail set at -5°, “3°, a nd 1° 
incidence. From these tests the trim angles of the horizontal tail 
were determined. (See table II.) 

Static tests . - In order to determine the static-stability 
derivatives C?,^, C n ^, and Cy^, the model was tested at i|r = ±5° 

through an angle -of -attack range of a ■= -2° to a = 23° for the flaps 
up (trim C L = O.33) and a = -2° to a = 18° for the flaps down 
(trim = O.76) for each of the slot configurations. Tests were also 

made at all selected test trim lift coefficients through an angle -of -yaw 
range of \|r = +20° to determine the variations of Cj, C n , and Cy 

with \|f for all slot configurations. 
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Tare teats were made for the 4-0 -per c ent - span alot configuration 
(flaps up and down) . The effect of the slots on the tares was assumed 
to he small; therefore, the tares for the 4-0-percent-span slot 
configuration were applied to all configurations. The Mach and Reynolds 
numbers for the tests were 0.17 and 1.01 x 10^, respectively. 

Rolling-flow tests . - Tests were conducted in the 6 -foot-diameter 
rolling-flow test section of the Langley stability tunnel, wherein 
rolling flight of the model was simulated by rotating the air stream. 

The model was mounted rigidly on a conventional support strut. Details 
of this test procedure are given in reference 1. 


All slot configurations were tested through the angle -of -attack 
range with the flaps up (trim Cl = 0.33) and with the flaps down 
(trim Cl = O. 76) at helix angles pb/2V of 0, ±0.0253, and ±0.0757 radian. 
The slopes of Cj, C n , and Cy plotted against pb/2V are the 


derivatives C 


V 


3^, and Cy^. The 40 -percent -span slot configuration 

(flaps up and down) was tested at the selected trim lift coefficients for 
the previously mentioned values of pb/2V from \|r = 0° to \|r = 20° 
to determine the variation of Ci , C n _, and Cy with \|r. The Mach 

r r Jr 


number and Reynolds number for the rolling -flow tests were 0.17 and 
1.01 x 10°, respectively. 


Yawing -flow tests . - Yawing-flow tests were conducted in the 
6- by 6-foot curved -flow test section of the Langley stability tunnel. In 
this section, curved flight is simulated approximately by directing 
the air in a curved path about a fixed model. 

All slot configurations were tested in curved flow through the 
angle -of -attack range with the flaps up (trim Cl = 0.33) and with the 
flaps down (trim Cl = O. 76) at values of rb/2V of 0, -0.039, -0.082, 
and -0.108. The slopes of C^, C n , and Cy plotted against rb/2V 

are the derivatives C j r , Cn r , and Cy r . 

The 40 -percent -span slot configuration (flaps up and down) was 
tested through an angle -of -yaw range of \|f = ±20° for values of rb/2V 
of 0, -0.039, -0.082, and -0.108 at the previously mentioned trim lift 
coefficients to determine the variation of C^ r , .C nr , and Cy r 

with \| 1 . The values presented herein are the average of the results 
at corresponding positive and negative angles of yaw. 

The 4-0 -percent -span slot configuration (flaps up) was tested at a 
trim lift coefficient of 0.33 through the angle-of-yaw range with the 
propeller off. 
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The yawing -flow tests were made at a Mach number of 0.13 
and a Reynolds number of 0.8 x 10°. 


Corrections 

Approximate jet-boundary corrections based on methods derived for 
unswept wings were applied to the angle of attack, longitudinal-force 
coefficient, and rolling -moment coefficient, and a blocking correction 
of 1.01 was applied to the dynamic pressure. 

Corrections for the effect of the support strut have been applied 
to C x , C L , Ci , C n , C Y , C V Cry and C Yr In rolling flow and 

curved flow, accurate tares were difficult to obtain, and as a result, 
the derivatives Ci , C np , C Yp , Ci r , C nr , and C Yr are not 

corrected for the effects of the support strut. 

The derivative Cj was corrected for the effective pitching 

velocity, which exists when the model is tested at an angle of yaw, by 
the following eqaution: 


C 1 = C-, ' cos i + f(A,Aj\|r) 

4 P 

where C i ’ is measured about the wind axis, and f(A,A , \|r), which is 
small as compared to C Zp ' cos t, is a function determined by use of the 

methods of reference k. Corresponding effective pitching corrections 
were not applied to the derivatives C np and C Yp . 

A correction was also applied to the derivative C Yr to account for 
the error caused by the cross -tunnel static -pressure gradient which is 
associated with curved flow. 


Experimental Results 

The experimental data are discussed briefly with reference to the 
effects of the slots and angle of yaw on the aerodynamic characteristics 
of the model, because the effects of these variables on the rotary 
derivatives have not been investigated extensively to date. The figures 
which present the results obtained in the present investigation are 
listed in table III. 
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The "basic lift and longitudinal -force data of figures 7 and 8 are 
generally in good agreement with larger scale tests of another model of 
the same airplane, as given in reference 7. 


The main effect of the slots is to extend the linear* range of those 
stability derivatives which are largely contributed by the wing to 
higher lift coefficients in a manner similar to the effect of slots on . 
the lift curve of a wing. (See figs i 7 to 10 and 12, 13, 15, and- l6.) 
One significant effect of the slots is on the damping in yaw C n ^ which 

increased as the slot span is increased. When the 80 -percent-span slots 
are used Ch is increased (over that of the unslotted configuration) 

by about 25 percent at = 0. (See fig. 15(a).) When the flapsare 
deflected (fig. 15(h)), the effects of the slots on the yawing-stability 
derivatives are not as great as when the flaps are retracted. 


A comparison of figures 11(a) to 11(f) shows that C n ^ varies to 


some extent through the yaw range, 
for any. given lift coefficient Cy 




Cy^ is approximately constant, and 
is approximately constant 


between i|r = -10° and \|r = 10°. As the angle of yaw is increased 


(fig.. 14(a)), 
constant, and 


C Y r 


tends to decrease, 


remains approximately 
decreases slightly. The variations are similar 


when the flaps are deflected. (See fig. 14(b).) 


The tests for the determination of the variation of Cy^, ( - : n r > 

and C i r with were made for negative values of rb/2V only. 

Results for positive values of rb/2V and positive \|r were obtained 
by assuming that the model was essentially symmetrical about the 
XY -plane and by utilizing the results for the corresponding opposite 
angles of yaw and rb/2V with regard for signs. This procedure amounted 
to averaging the derivatives for corresponding positive and negative 
angles of yaw. 

In general, as the angle of yaw is increased from 0° to 20°, the 
damping in yaw C nr is increased by about 15 percent, and Cy r 

and are decreased slightly (fig. 17(a)). Deflection of the flaps 

or removal of the propeller does not appreciably change the variations 
of the yawing -stability derivatives, with angle of yaw. (See fig. 17). 

Because C nr is largely dependent on the size of the vertical tail 
it is approximately true throughout the yaw range that 

2 -K 
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Substitution of the proper values of the span, 'C n ^, and. the tail 
length ly indicates that the trend of the variation of C nr with 
‘Shown in figure 17 is reasonable. 


The relative constancy of C Yp , C lp , C np , C Yr , C lr , and C nr 

with angle of yaw as indicated by figures 14 and 17 and the linearity of 
the curves of C Y , C^, and C n plotted against \|r for angles of yaw 

up to approximately 10° (fig. 11) were factors which indicated that 
nonlinearities were not of first-order importance for this airplane in 
the calculation of motions involving reasonably small variations in \|r. 
Consequently, most of the motion calculations neglect the effect of \|r 
on the stability derivatives. 


The results of unpublished tests of swept wings at Reynolds numbers 
to 8.0 x 10° in the Langley 19-foot pressure tunnel indicate that the 
linear part of the curve. of Cj plotted against is increased by 

an increase in Reynolds number. The curves of Cj ^ against given 

herein agree well with those obtained in flight (references 5 and 6) 


and in tests of a •jj^r-scale model (reference 7) except at lift coefficients 
above Cp = 0.8 where the magnitude of the present test values decrease 


for the slotted -wing configurations. The linear parts of the curves 
of the rolling- and yawing-stability derivatives are believed to be 
extended similarly to higher lift coefficients if the Reynolds number 
is increased. 


MOTION CALCULATION METHODS 


. The lateral disturbed motions of the test airplane were calculated 
from the aerodynamic data obtained from the tests described in the 
section entitled "Wind-Tunnel Tests." The mass and dimensional 
characteristics of the airplane are given in table I, and a tabulation” 
of the flight conditions for which lateral disturbed motions were 
calculated for comparison with the flight motions are given in table 17. 
Most of the calculations involved dynamic derivatives which were constant 
for a given lift coefficient as is usually employed in the theory of 
small disturbances used in lateral -stability calculations. 

Solutions of the lateral equations of motion, given in the appendix, 
were obtained. for unit step disturbances in roll and yaw by the method 
described in reference 8. The aileron and rudder deflections during the 
flight motion under investigation were then approximated by a series 
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of step functions usually at ^--second intervals. Tlie motion arising 
from the control movements was then calculated by 

(&m-n) rudder (Bn) yaw + ( &m-n) aileron^ ^ roll 


m 


P n > Pn> or r n = 


n=0 


where is the value of the unit solution caused "by AC n = 1 

or AC 1 =1 at a- time T = kn, and &m-n is the fraction of the unit 
disturbance applied by the rudder or aileron at a time T = k(m - n) . 

The rudder and aileron effectiveness were obtained from reference 7- 
This procedure is essentially an approximate evaluation of Duhamel's 
integral and was considered sufficiently accurate for these calculations. 
Reference 8 gives a more exact graphical evaluation of this integral. 

The yawing moment caused by aileron deflection and the rolling moment 
caused by rudder deflection were not considered of enough significance 
to warrant their inclusion in these calculations. In some cases, 
however, these factors mhy be of greater significance. 

In a few cases, unit solutions to the equations of motion were 
obtained by a Laplace transform procedure which has been adapted for 
use with automatic digital computers. The results were, of course, 
identical with those presented in this paper. 

Calculations of the lateral motions for a few cases employing a 
nonlinear variation of rolling-moment and yawing-moment coefficients 
with angle of sideslip and a variation of Cn r with angle of sideslip 

were carried out by use of the Kutta three-eighths rule for solving the 
lateral equations of motion. (See reference 9- ) All lateral -mot ion 
calculations were made on an automatic digital computing machine. 


CALCULATED LATERAL MOTIONS 
General 


The flight records corresponding to the motions calculated for this 
paper showed that the motions resulting from right and left control 
movements were not exactly of opposite magnitude. This result indicates 
that there was some asymmetry in the characteristics of the airplane, 
although the wind-tunnel test results indicated no marked asymmetry in 
the model characteristics. Part of the differences in the flight 
motions to the right and left are undoubtedly due to variation in the 
thrust conditions. Although the flight tests were made under approxi- 
mately zero thrust, no convenient means of obtaining this condition was 
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available, and thus this requirement was left to the judgement of the 
pilot. Some small part of the difference in the motions to the left 
and right should be attributable to instrument error as the film record 
frequently contained considerable hash which, of course, tended to 
obscure the actual motion record. The differences between the motions 
to the right and left were resolved in the present paper by presenting 
both records with the sign .of one reversed so that the motions were 
superimposed. The shaded areas in the lateral-motion plots (figs. 18 
to 23 and. 25 to 35) represent the difference between the motions to the 
right and left. The flight motions corresponding to the calculated 
motions were obtained from references 5 and 6 and related tests. The 
figures which present the results of the lateral -motion calculations are 
listed in table III. 


Lateral Oscillations 


Lateral oscillations were initiated, by abruptly deflecting the rudder 
of the airplane and returning it to neutral equally as rapid a moment 
later. Some aileron waggle caused by the floating tendency of the 
ailerons occurred and was accounted for in the calculations by the 
procedure given in the calculation methods. 


The calculated and flight motions are generally in good agreement 
for all motions calculated (figs. 18 to 23) except for the condition 
at = 0.977* (See fig. 23.) A comparison of the flight and 

calculated periods and time to damp to half- amplitude (fig. 2h) also 
indicate fairly good agreement except for the motions at lift coefficients 
of 0.977 and 1.169. 

The values of used for the calculation of the lateral motions 

were obtained from the curves of plotted against ^ (fig. 11) 

wherever possible rather than the curves of Ci obtained from tests 

♦ 

made at = +5° (figs. 9 and 10) because this procedure is believed 
to be more accurate. Although the difference between the two methods 
is generally small such is not always true, and in some cases the 
calculated rate of damping was found to be appreciably affected by the 
difference in Ci n . 


At the higher lift coefficients, the experimental stability 
derivatives deviate appreciably from their initial trends. 'This 
tendency previously has been referred to as a Reynolds number effect 
and is probably the cause for the lack of agreement at high lift 
coefficients between the flight and calculated results. References 5 
. and 6 which present flight tests for this airplane show no reduction 
in Cj up to the maximum test lift coefficient. The rotary derivatives 
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of the airplane presumably hehave similarly. Evidence of the deviation 
of the experimental stability derivatives from their true variation with 
lift coefficient is observed in figure 20 which presents results for 
a lateral oscillation occurring at = 0.759- The calculated result 

indicates an excessive response in sideslip, yaw, and roll. 

A few additional calculations of the period and time to damp, made 
for other high lift coefficients in flight, indicate increasingly poor 
agreement between the flight and calculated times to damp to one -half 
amplitude with increasing lift coefficient. The relatively good 
agreement between the periods of the flight and calculated motions for 
all lift coefficients, however, indicates that the experimental values 
of C n „ are fairly close to the correct values. The period of the 
motion is primarily a function of the directional- stability para- 
meter C n p. (See reference 10.) An extrapolation of the curves of the 

derivatives plotted against lift coefficient - which amounted to selecting 
the value of the derivative just previous to the break in the curves 
occuring at maximum lift coefficient - was employed for one case' 
at C L = 1.17 "but failed to yield a satisfactory result. A linear 
extrapolation of the curves in the region preceding the departure from 
the theoretical or linear trend is expected to be more satisfactory. 


Eudder Kicks 

Eudder kicks were initiated by abruptly deflecting the rudder of 
the airplane and holding this deflection as the airplane responded. 

The records of the flight motions were short, usually covering 5 seconds. 
As in the case of the lateral oscillations discussed previously, Borne 
slight aileron waggle that occurred was accounted for in the calculated 
motions . The agreement between the flight and calculated motions is, 
in general, quite good (figs. 25 to 3l) • Agreement is best at the 
lower lift coefficients; however, at a lift coefficient of 0.919 (fig. 28) 
good agreement is shown between flight and calculated angles of sideslip 
and rolling velocities. 

The flight tests reported in reference 5 showed that at low lift 
coefficients the airplane rolled in response to a rudder kick as if 
the airplane had a decided negative dihedral effect. Figure 31 
indicates that this peculiar response was largely caused by the slight 
aileron waggle which occurred during the motion. In general, for all 
configurations investigated, it was found necessary to account for any 
slight aileron movements in order to satisfactorily predict the lateral 
motions of the airplane. 

The motions at C L = 0.79^ (fig- 30) give evidence of departure of 
the experimental derivatives from the true variation with lift coefficient 
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This result again is the Reynolds number effect previously discussed. 
The flight and calculated periods of the lateral oscillation caused by 
the rudder kicks are in very good agreement which indicates again that 
the values of C n used in the- calculations were reasonably accurate . 


Aileron Rolls 

Aileron rolls were initiated by abruptly deflecting the aileron 
control of the airplane and holding the deflected position as the 
airplane responded. The rudder control was held as near neutral as 
possible. The records of motion were necessarily short because of the 
large angles of bank assumed by the airplane after a short period of 
time. 


Comparison of the flight and calculated aileron rolls indicates 
fairly good agreement (figs. 32 to 35) except for the high-lift-coefficient 
condition without nose slots. The stability derivatives for the 
unslotted configuration show a departure at fairly low lift coefficient 
from the initial trend of the variation of the derivatives with lift 
coefficient. As previously mentioned, this result is not obtained at 
the flight Reynolds numbers. The calculations for this high lift 
coefficient were carried out with two sets of stability derivatives. 

One set was obtained by a linear extrapolation of the curves of the 
derivatives plotted against lift coefficient; and the second set by 
selecting for the value of the derivative that value Just previous to 
the final break in the curves of the derivatives plotted against lift 
coefficient. The second procedure was necessary because the flight 
lift coefficient was greater than the maximum experimental lift 
coefficient. The linear extrapolation produces the better result, but 
it can be seen that neither set of derivatives was very satisfactory, 
although it would appear logical to believe, in view of the Reynolds 
.number effects indicated previously, that the linear extrapolation of 
the derivatives would have given fairly good results. The aileron 
roll for the 80 -percent -span slot configuration with flaps down 
at C L = I.169 (fig. 35) shows fair agreement between the flight and 
calculated motions so that a beneficial effect of the slots in 
maintaining an unseparated flow over the wing to a fairly high lift 
coefficient is indicated. 

The effects of elasticity of the wing were not considered in the 
calculations of the aileron rolls because of lack of knowledge of the 
flexibility of the wing. Rough estimates of the effect of elasticity 
on the rolling velocity for the results of figure 33 indicate that the 
discrepancies shown between flight and calculated results would be 
reduced by about 50 percent or more . 
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A comparison has been made of the maximum rolling velocities 
obtained by calculation for conditions duplicating those in flight, by 
flight tests, and by calculations for a coordinated maneuver in. which 
the yawing velocity and angle of sideslip are maintained at zero. For 
this last case, the rolling velocity becomes 


P 


-AC 


I i V 
b 


The results are presented in the following table: 


Figure 

Maximum rolling velocities from figures 


% 

2V 

b 

Calculated 

Flight 

(left) 

Flight 

(right) 

32 

0.62 

0.44 

— 

2.67 

33 

• 77 

.60 

0.53 

1.11 

3^ 

• 78 

— 

.67 

1.13 - 

"35 

.34 

.33 

.30 

.68 


The comparison indicates that a marked difference inay be obtained in 
the maximum rolling velocity by eliminating the degrees of freedom of 
sideslip and yaw. This last method of calculation is used frequently 
in comparing the relative merits of various forms of ailerons. 


Nonlinear Calculations 

Calculations in which the curves of C^ and C n against t were 

represented not by a single slope but by a series of tangents to the 
curves, and in which the variation of C n ^ with \|r was included were 

made for two motions to illustrate the effect of these nonlinearities. 
The results are given in figures and 37 along with the flight motion 
and the calculations made for linear slopes and constant damp ing 
derivatives. The effect of the aileron waggle was small and was taken 
to be the same as for the linear calculation. These figures ind icate 
that for the angles of sideslip encountered in these motions little 
is to be gained by going to the additional effort required to make 
the calculations for the nonlinear case. The nonlinear calculations 
required approximately ten times as long to complete as the linear 
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calculations. Motions at large angles of sideslip or. motions for a 
configuration having’ more pronounced nonlinearities than for the 
airplane co ns idered herein would undoubtedly show the nonlinearities to 
he of greater significance. 


SUMMARY OF RESULTS 


The investigation indicated that the lateral disturbed motions of 
the airplane considered herein may be generally calculated with good 
accuracy by the use of experimental-scale -model data up to lift 
coefficients where any appreciable effects of the difference in Reynolds 
number between the scale and flight conditions begin to be apparent. 

This effect is usually evidenced by a departure of the scale results 
from a gradual variation with lift coefficient and, for the case 
investigated, corresponds to a lift coefficient of 0.8 for the slotted- 
wing configurations. At higher lift coefficients the accuracy in 
calculating flight motions is progressively poorer, especially in the 
time to damp to half- amplitude where errors as much as 100 percent or 
more may be encountered. For the airplane of this investigation, at 
these high lift coefficients, the calculated time to' damp to half 
amplitude was generally higher than the flight value. In most cases 
the period of the motion could be calculated to within 5 percent of that 
obtained in the flight motion. In order to obtain good accuracy in 
predicting flight motions proper consideration of all control movements 
must be made. 

For most cases, the nonlinear variation of the directional- stability 
parameter, the effective dihedral parameter, and the variation of the 
damping in yaw with angle of yaw could be neglected. If, however, the 
flight motions reach a sufficiently large amplitude for a long period 
of time, as in an aileron roll for example, greater accuracy in the 
calculated motions should result if the nonlinearities in the variation 
of yawing-moment coefficient, rolling -moment coefficient, and damping 
in yaw with angle of yaw are included in the calculations. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Ya., November 1, 19^9 
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APPENDIX 

< 

LATERAL EQUATIONS OF MOTION 


The lateral equations of motions used for most of the calculations 
of this paper consider the fluctuation of the aerodynamic forces to tie 
directly proportional to the component angular and translational 
velocities associated with the axis system as is customary in classical 
stability theory. The modification to these equations to account for 
the nonlinear effects of sideslip are discussed in the text. The 
equations used herein contain the necessary product of inertia and 
control terms and are given as follows: 


Roll 


Cz 5r S r + C z & a + C Zp 3 + |c Zp Dcp + = 2 J 
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—I CO^r, 
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^0 I 0 1 
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cos q sin q D 2 ^ ^ 
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where 


D= i 


TV 


8 = ■ 
C l6 r 

° ls a 

Cn 8 r 

Cn &a 

° Y & 

u r 

°Y 8 a 

<P 

d = 

W 

g 

k Xo 

k Z D 

7 

V 


rate of change of rolling-moment coefficient with rudder 
-deflection 

rate of change of rolling -moment coefficient with aileron 
deflection 

rate of change of yawing -moment coefficient with rudder 
deflection 

rate of change of yawing -moment coefficient with aileron 
deflection 

rate of change of lateral-force coefficient with rudder 
deflection 

rate of change of lateral-force coefficient with aileron 
deflection 


W 

pStg 


angle of "bank, radians 

weight of airplane 
acceleration due to gravity 

radius of gyration about principal longitudinal axis 


radius of gyration about principal normal axis 


flight path angle j positive for climb 


(/ 
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inclination of principal longitudinal axis of inertia with 
respect to flight path; positive when above flight path 
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TABLE I.- PERTINENT AIRPLANE DIMENSIONS AND CHARACTERISTICS 
Mass characteristics: 

Normal gross weight, lb_ 8,700 

Moments of inertia, ft-lb-sec2 : 

Ix 0 7,654 

• Iy 0 14,088 

*Z 0 20,159 

Principal axis (relative to thrust line at nose) 0.45° up 

Center of gravity: 

Location on M. A. C. , percent M. A. C. ...... . . . . . 21.8 

Relative to thrust line, percent M.A.C. below 14.8 

Wing: 

Spar., ft 33.6 

Area, sq ft 250 

Airfoil section (normal to L.E. ) : 

Root . Modified 66,2x-ll6 (a=0.6) 


M.A.C. , ft 7-79 

Leading edge of M.A.C. (ft behind L.E. root chord) 3-27 

Aspect ratio 4.51 

Taper ratio 1.84:1.00 

Dihedral, deg 0 

Sweepback (at quarter -chord line), deg . 35 

Total area, plain, sealed wing flaps, sq ft 12.52 

Ailerons : 

Span (along hinge line, each), ft 8.75 

Area (rearward of hinge center line, each), sq ft 6.51 

Horizontal Tail: 

Total area, sq ft 46.53 

Vertical tail: 

' Fin area (above horizontal tail), sq ft 13-47 

Total rudder area, sq ft 10.26 

Ventral fin area, sq ft: 

. Large ventral 17-10 

Small ventral (approx.) 8.50 
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TABLE II.- TRIM ANGLES OF HORIZONTAL TAIL 
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TABLE III . - INDEX OF FIGURES 

Wind-tunnel results: 

Static -stability characteristics: -Figures 

Longitudinal characteristics 7 and 8 

Lateral derivatives 9 and 10 

Variations of C 1 , C n , and Cy with ijr 11 

Boiling derivatives : 

Derivatives at i|r = 0° . 12 and 13 

Variations of derivatives with ’ 14 

Yawing derivatives: 

Derivatives at \|/ = 0° .15 and 16 

Variation of derivatives with i|f .' . . . .17 

Calculated and flight motions: 

Abrupt deflection and release of rudder 18 to 23 

Summary of period and damping values 24 

Eudder kicks 25 to 31 

Aileron rolls 32 to 35 

Nonlinear aerodynamic effects of sideslip 36 and 37 
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TABLE IV. - FLIGHT CONDITIONS FOR WHICH MOTIONS WERE CALCULATED 


Figure 

Slots 

(percent span) 

Flaps 

(d.eg) 

Ventral 

fin 

Vc x 

(mph) 

C L 

Type of test 

18 

40 

0 




198 

0.334 



19 

4o 

0 




156 

• 551 



20 

40 

0 


>Large < 


136 

• 759 


>Oscillation 

21 

40 

40 




160 

.524 



22 

40 

40 




128 

.801 



23 

80 

0 


Small 


120 

.977 



25- 

40 

0 

"S 



198 

.341 


. 

26 

40 

0 




160 

.556 



27 

40 

0 




133 

.764 



28 

40 

0 




120 

•919 

; 

>RuddLer kick 

. 29 

40 

40 


^ Large < 


157 

.537 



30 

40 

40 




130 

.794 



31 

80 

0 




228 

.278 



32 

0 

0 




119 

.983 

> 


33 

40 

0 




150 

.598 


>Aileron roll 

34 

40 

40 

I 



146 

.598 



35 

8o 

'40 j 

J 



110 

1.169 

-> 
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Figure 1.- System of stability axes. Arrows indicate positive forces, 

moments, and angular displacements. 
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Langley stability tunnel. 
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Figure 5.- Close-up of UO-percent leading-edge slots on Q-scale model. 
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Figure 7*- Concluded. 
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Figure 10. - Variation of lateral static- stability derivatives with lift 
coefficient for two slot configurations. Trim C L = 0.33; t = 0°; 

8f = 0°; R = 1.01 X 10 6 . 
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(a) 0-percent -span slots; = 0°. . 

Figure 11.- Variation of Cy> C n , and with angle of yaw. 

ventral fin on; propeller on; R = 1.01 X 10°. 
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Cb) O-percent- span slots; 8p = U0°. 
Figure 11.- Continued. 
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(e) 80- percent -span slots; B-f = 0°. 
Figure 11.- Continued. 
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(g) J+O-percent-span slots; Sf = 0°; propeller off. 
Figure 11.- Continued. 
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Lift coefficient , 


(a) Trim C L = 0.33; B f = 0°. 

Figure 12.- Variation of rolling- stability derivatives with lift 
coefficient for three slot configurations, = 0°; 

R = 1.01 x 10 6 . 
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Figure 13.- Variation of rolling- stability derivatives with lift 
coefficient for two slot configurations, t = 0°j = 0°; 

R = 1.01 x 10 6 . 
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• (a) 5 f = 0°. 

Figure lb.- Variation of rolling- stability derivatives with angle of 
yaw. Large ventral fin on; propeller on; ^-percent- span slots; 

R = 1.01 x 10 6 . 
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Figure 15* - Concluded. 
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Figure 17.- Variation of yawing- stability derivatives with angle of yaw 
Large ventral fin on; 40-percent-span slots; R = 0.8 x 10°. 
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Time, see 


Figure 20.- Comparison of flight and calculated lateral motions 
resulting from an abrupt deflection and release of the rudder, 
to-percent-span slots; flaps up; Cl = 0.759; Vc = 136 miles per hour 
engine idling. 
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Figure 24.- Comparison of flight and calculated values of the period and time to damp to 

half amplitude. 
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Figure 26.- Comparison of flight and calculated lateral motions 
resulting from a left rudder kick. 1+O-percent-span slots; 
flaps up; = 0.556; V c = 160 miles per hour; engine idling. 
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Figure 28.- Comparison of flight and calculated lateral motions 
resulting from a right rudder kick. 1+0-percent- span slots; 
flaps up; C L = 0.919; V c = 120 miles per hour; engine idling. 
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Figure 29.- Comparison of flight and calculated lateral motions 
resulting from a left rudder kick. 40-percent-span slots; 
flaps down; C L = 0.537; V c = 157 miles per hour; engine idling. 
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Figure 30.- Comparison of flight and calculated lateral motions 
resulting from a left rudder kick. 40-percent-span slots; 
flaps down; C^ = 0.79^; V c * 130 miles per hour; engine idling. 
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Figure 32.- Comparison of flight and calculated lateral motions 
resulting from an abrupt left aileron deflection. O-percent- 
span slots; flaps up; = 0.983; V c = 119 miles per hour; 

engine idling. 
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Figure 33-- Comparison of flight and calculated lateral motions 
re siilting from an abrupt left aileron deflection. 1+0 -percent 
span slots; flaps up; C-^ = 0 . 598 ; V c = 150 miles per hour; 

engine idling. 
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Figure 3^*- Comparison of flight and calculated lateral motions 
resulting from an abrupt right aileron deflection. 40-percent- 
span slots; flaps down; C L = 0.598; V c = 1^6 miles per hour; 
engine idling. 
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Figure 35*- Comparison 
resulting from an a 1 
span slots; flaps d 

engine idling. 
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Figure 36 .- Effect of nonlinear variation of C n , and with 

angle of yaw on the calculated lateral motions resulting from a 
left rudder kick. 40-percent-span slots; flaps up; = 0.764; 

V c = 133 miles per hour; engine idling. 
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Figure 37* - Effect of nonlinear variation of Cj, C n , and C n with 

angle of yaw on the calculated lateral motions resulting from an 
abrupt left aileron deflection. 1*0 -percent -span slots; flaps up; 
C L = 0.598; V c = 150 miles per hour; engine idling. 
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